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We discuss the production of hadronic resonances in very peripheral heavy ion collisions, where the 
ions collide with impact parameter larger than twice the nuclear radius and remain intact after the 
collision. We compare the resonance production through two-photon and double Pomeron exchange, 
showing that when we impose the condition for a peripheral interaction the 77 process dominates 
over the Pomeron interaction, due to the short range propagation of this last one. We also discuss 
the observation of light resonances through the subprocess 77 i? ^ 77, which is a clean signal 
for glueball candidates as well as one way to check the existence of a possible scalar a meson. 



I. INTRODUCTION 

Collisions at rclativistic heavy ion colliders like the 
Relativistic Heavy Ion Collider RHIC/Brookhaven and 
the Large Hadron Collider LHC/CERN (operating in its 
heavy ion mode) are mainly devoted to the search of 
the Quark Gluon Plasma. In addition to this important 
feature of heavy-ion colliders, we will also have ultra- 
peripheral collisions with impact parameter b > 2Ra, 
where Ra is the nuclear radius, and where the ions re- 
main intact after the collision. 

These interactions will be mostly of electromagnetic 
origin: two-photon (77) or photonuclear processes (7^). 
Due to the very strong photon field of each charge Z 
accelerated ion, the photon luminosity will be quite high. 
In the case of RHIC final states produced in the two- 
photon process with an invariant mass up to a few GeV 
will appear at large rates. Above this scale the photon 
luminosity drops very fast. At LHC a final state with a 
mass almost two orders of magnitude larger can still be 
produced at reasonable rates. The variety of processes 
that can be studied in heavy ion peripheral collisions have 
been extensively reviewed recently [ p . 

The fact that hadronic resonances (i?) could be pro- 
duced at large rates in peripheral heavy ion collisions was 
already discussed many years ago Perhaps this may 
be one of the most interesting studies to be performed 
at RHIC, because the machine will serve as a factory of 
light hadrons in 77 and jA reactions. Vector resonances 
will appear at huge rates in photonuclear reactions |^ . as 
well as scalar and pseudoscalar resonances in two-photon 
processes Q. We will particularly focus our attention 
on the production of scalar and pseudoscalar resonances 
through the 77 process. 

Two-photon physics at e+e^ colliders provided for a 



long time a lot of information on hadronic resonances 
The two-photon process is very important because 
it involves the electromagnetic coupling of the resonance, 
and its knowledge with high precision is very useful, for 
instance, to unravel the possible amount of mixing in 
some glueball candidates ||^ , complementing the informa- 
tion obtained through the observation of hadronic decays. 
Another interesting study is the possible production of a 
light scalar meson (a) whose existence has been discussed 
for several years |^]. We stress again that the advantage 
of relativistic heavy ion colliders is that the photon lu- 
minosity for two-photon physics is orders of magnitude 
larger than the one at available e~^e~ machines. 

We will discuss the production of light hadronic reso- 
nances in ultra-peripheral heavy ion collisions. We will 
show that the process 77 ^ i? 77. see Fig.(|), can 
be observed for many resonances above or at the same 
level of the background. The main background is the 
continuum reaction 77 — > 77, this one will be discussed 
as well as some other background contributions. Dou- 
ble Pomeron exchange may also compete with the 77 
physics, we will point out that this contribution is not 
important for very heavy ions. 




* Talk at the XXII Brazilian National Meeting on Particles 
and Fields, Sao Lourengo, Minas Gerais, 2001 



FIG. 1. Diagram for 77 fusion in a peripheral heavy-ion 
collision. The blob represents a continuum or resonant process 
leading to a two-photon final state. 
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The distribution of this review is the following: In Sec- 
tion II we present the distribution functions for photons 
and Pomerons in the ion, and discuss some of the approx- 
imations to obtain realistic cross sections. We compare 
77 processes with the ones initiated by double Pomcron 
exchange. In Section III we discuss the 77 ^ i? ^ 77 
reaction and its background in the case of some glueball 
candidates and in the case of a possible scalar a meson. 
Section IV contains our conclusions. 



II. TWO-PHOTON AND DOUBLE POMERON 
EXCHANGE PROCESSES 

A. Distribution functions 

The photon distribution in the nucleus can be de- 
scribed using the equivalent-photon or Weizsacker- 
Williams approximation in the impact parameter space. 
Denoting by F{x)dx the number of photons carrying a 
fraction between x and a; -I- dec of the total momentum 
of a nucleus of charge Ze^ we can define the two-photon 
luminosity through 

^= f^F{x)F{T/xl (1) 

where r = s/s, s is the square of the center of mass 
(c.m.s.) system energy of the two photons and s of the 
ion-ion system. The total cross section of the process 
AA -> ylyl77 is 

, , f dL , , 
'^(•^) = J dr^^{§), (2) 

where a{s) is the cross-section of the subprocess 77 — > X. 

There remains only to determine F{x). In the liter- 
ature there are several approaches for doing so, and we 
choose the conservative and more realistic photon distri- 
bution of Ref. Q . Cahn and Jackson Q , using a prescrip- 
tion proposed by Baur , obtained a photon distribution 
which is not factorizable. However, they were able to give 
a fit for the differential luminosity which is quite useful 
in practical calculations: 

where z — 2MR^, M is the nucleus mass, R its radius 
and ^(z) is given by 

C(^) = ^A,e-^•^ (4) 
1=1 

which is a fit resulting from the numerical integration 
of the photon distribution, accurate to 2% or better for 
0.05 < z < 5.0, and where Ai = 1.909, A2 = 12.35, 



A3 = 46.28, hi = 2.566, &2 = 4.948, and 63 = 15.21. For 
z < 0.05 we use the expression (see Ref. g) 

The condition for realistic peripheral collisions (&min > 
Ri + R2) is present in the photon distributions showed 
above. 

The processes that we shall discuss can also be inter- 
mediated by the diffractive subprocess PP X, where 
P is the Pomeron. 

In the case where the intermediary particles exchanged 
in the nucleus-nucleus collisions are Pomerons instead of 
photons, we can follow closely the work of Miiller and 
Schramm ]lO[ | and make a generalization of the equiv- 
alent photon approximation method to this new situa- 
tion. So the cross section for particle production via two 
Pomerons exchange can be written as 

O-AA = j dxidX2fp{xi)fp{x2)(Jpp{spp), (6) 

where fp{x) is the distribution function that describe 
the probability for finding a Pomeron in the nucleus with 
energy fraction x and <Tpp{spp) is the subprocess cross 
section with energy squared s pp . In the case of inclusive 
particle production we use the form given by Donnachie 
and Landshoff ||ll| 

fp{^) = ^ dt\f3Ap{t)\^ \Dp{t-s')\\ (7) 

where Dp{t;s') is the Pomeron propagator ||l^ 

sm(i7rap(t)) V 2 / 

with s the total squared cm. energy. The Regge trajec- 
tory obeyed by the Pomeron is ap{t) = 1+e + a'pt, where 
e — 0.085, a'p — 0.25 GeV~^ and i is a small exchanged 
four-momentum square, t = fc^ << 1, so the Pomeron 
behaves like a spin-one boson. The term in the denom- 
inator of the Pomeron propagator, [sm('^Trap(t))]~^ , is 
the signature factor that express the different properties 
of the Pomeron under C and P conjugation. At very high 
cm. energy this factor falls very rapidly with k^, whose 
exponential slope is given by a'p ln(s/m^), m is the pro- 
ton mass, and it is possible to neglect this dependence, 

sm^TT{l + e - a'pk^) fa cos{^'Ke) fa 1. (9) 

If we define the Pomeron range parameter vq as 

rl^a'p\n{s/m^), (10) 
the Pomcron propagator can be written as 

\Dp{t = -k^; s)| = {s/mye~'-''^\ (11) 
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Since we are interested in the spatial distribution of the 
virtual quanta in the nuclear rest frame we are using 

t = -k2. 

The nucleus-Pomeron coupling has the form 



l3Ap{t)^3ApQFAi-t), 



(12) 



where /3o = 1-8 GeV~^ is the Pomeron-quark coupling, A 
is the atomic number of the colliding nucleus, and t) 
is the nuclear form factor for which is usually assumed a 
Gaussian expression (see, e.g., Drees et al. in 



FA{-t) = e'/^Qo^ 



(13) 



where Qo = 60 MeV. 

Performing the t integration of the distribution func- 
tion in Eq.(|7|) we obtain 



(3v4/3o)- 
(27r)2a; 
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The total cross section for a inclusive particle produc- 
tion is obtained with the above distribution and also with 
the expression for the subprocess PP — X as prescribed 
in Eq.d). 



taking into account the impact parameter dependence, 
b. It is straightforward to verify that in the collision of 
two identical nuclei the total cross section of Eq.(|^) is 
modified to 



^pp- 



Q'2 



d26 



27r 



b /2 PP 
^AA^ 



(14) 



2rg . The total cross section for 



where (Q')"' = iQoY 
inclusive processes is obtained after integration of Eq. (|14| ) 
with the condition bmm > '2Ra in the case of identical 
ions. 

Another way of to exclude events due to inelastic cen- 
tral collisions is through the introduction of an absortion 
factor computed in the Glauber approximation . This 
factor modifies the cross section in the following form 



da 



AA 



da^^- 



d^b 



d^b 

^ PP^R 

d^b 



exp 



'A^bao 



dQ^ 



i2ny 



JQb 



exp 



^A'bao ^ e 
47r 



(15) 



where CTq is the nucleon-nucleon total cross section, whose 
value for the different energy domains that we shall con- 
sider is obtained directly from the fit of Rcf. [n5| 



do = Xs"-|-Yis-''i +Y2S 



-V2 



(16) 



B. Is double Pomeron exchange a background for 77 
processes? 

The double Pomeron exchange producing a final state 
X have matrix elements with the same angular struc- 
ture as the 77 one. This is easy to observe within the 
Donnachie and Landshoff model ||l^. For instance, in 
this model to compute the cross section of the subpro- 
cess PP ^ i? it is assumed that the Pomeron couples 
to the quarks of the resonance (R) like a isoscalar pho- 
ton . This means that the subprocess cross section of 
PP R can be obtained from suitable modifications on 
the cross section for 77 i?, and it will differ only by 
the coupling constant and a form factor describing the 
phenomenological Pomeron-quark coupling. Therefore it 
is natural to ask if the Pomeron process is a background 
to the two-photon process, and when it has to be added 
to the calculation of a specific process. 

In general the double Pomeron exchange is not a back- 
ground for the two-photon process and this is easy to 
understand. The Pomeron contrarily to the photon does 
not propagate at large distances, actually its propagator 
has a range parameter tq defined in Eq.(^o|). When we 
impose the condition for ultra peripheral collisions {i.e. 
the nuclei do not physically collide) the cross section di- 
minishes considerably. 

In Eq.d) the cases where the two nuclei overlap are not 
excluded. To enforce the realistic condition of a periph- 
eral collision it is necessary to perform the calculation 



with X ^ 18.256, Yi = 60.19, Y2 = 33.43, e = 0.34, Tyi = 
0.34^2 = 0.55, Fa{Q^) = e-'S'/^Q^ and we exemplified 



pp- 

'AA 



R 



Eq.(|l5[) for the case of resonance production, i.e., a 
is the total cross section for the resonance production to 
be discussed in the sequence. The integration in Eq.(|T5|) 
is over all impact parameter space 

We compared the rates for double Pomeron exchange 
and two-photon production of several final states like res- 
onances, a pair of pions and a hadron cluster of invariant 
mass Mx- The details of the calculations can be found 
in Ref. and here we will describe part of the results 
of that work. 

In Table I we compare the cross sections for resonance 
production through the processes 77, PP R. 
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TABLE I. Cross sections for resonance production through 
photon-photon (77) and double-Pomeron (PP) processes. 
For RHIC, = 200 GeV/nucleon, we considered ^^*U ion 
and for LHC, y/s = 6300 GeV/nucleon, the nucleus is ^""^Pb. 
The cross sections are in mbarn. Rates computed with the 
geometrical cut 6 > 2Ra 
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Table II shows the ratios of cross sections for diffrac- 
tive resonance production calculated with the Glauber 
absorption factor to the one with the geometrical cut. 
The exclusion of central collisions through the Glauber 
absorption factor is stronger than the one with the ge- 
ometrical cut. Table III shows the 7r° production for 
different ions and at different energies. We observe that 
the rates for double Pomeron exchange becomes closer to 
the two-photon rates when we go to lighter ions. 



Meson 



V 
V' 

Vc 
V'c 

Vb 



3.54 X 10'"' 
3.58 X 10"^ 

3.46 X 10"^ 

3.47 X 10"^ 
3.61 X 10"^ 
3.5 X 10"^ 



1.5 X 10"^ 
1.47 X 10"^ 
1.5 X 10"^ 
1.32 X 10"^ 
1.5 X 10~^ 
1.45 X 10"^ 



TABLE II. Ratios of cross sections for diffractivc reso- 
nance production calculated with the Glauber absorption fac- 
tor to the one with the geometrical cut in the collision of ^-^^^U 
for energies available at RHIC (y^ — 200 GeV/nucleon), and 
collisions of ^'^^Pb for energies available at LHC (-^s — 6.300 
GeV/nucleon). 

The production of a cluster of particles through double 
Pomeron exchange is also dominated by the 77 process 
(see Fig.(l) of Ref. Q). 
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Si (A=28) 


100 


0.22 X 10"^ 


0.12 X 10" 


-3 


0.39 X 10"^ 



TABLE III. Cross section for vr" production for differ- 
ent ions and at different energies. The energies are in 
GeV/nucleon and the cross sections in mbarn. a^^~'^ is the 



cross section computed with the geometrical cut and a 
the one with the absorption factor. 



In general we may say that for very heavy ions the dou- 
ble Pomeron exchange gives cross sections one order of 
magnitude (or more) below the two-photon process. This 
changes when we collide lighter ions. As we go down from 
large charge ions to smaller ones the effect of Pomeron 
physics increases and it dominates the electromagnetic 
physics, as happens in the proton case. The fact that 
the Pomeron has a short range parameter is not the only 
fact that counts when we analyze each specific process. 
It must be also remembered that the Pomeron couples to 
light and heavy quarks differently. Apart from kinemat- 
ics the differences in the rates of resonance production 
by photons and Pomerons increases between those res- 
onances formed by light and heavy quarks, as seen in 
Table I. 



III. THE REACTION 7 + 7^7 + 7 

A. The continuum process 

We are particularly interested in the photon-photon 
scattering because it can be a very clean signal for 
hadronic resonances like glueballs and the a meson. On 
the other hand this scattering is important by itself, and 
could probably be directly measured by the first time at 
RHIC, as predicted in Ref. 

The subprocess 77 77 up to energies of a few GeV is 
dominated by the continuous fermion box diagram, and 
is a background for the resonant 77 — > i? — > 77 process. 
It was first calculated exactly by Karplus and Neuman 
[]l7| and Dc ToUis ||lj]. There are sixteen helicity ampli- 
tudes for the process and, due to symmetry properties, 
the number of independent amplitudes will be only five, 

that may be chosen to be , | , 

or — denotes the 



and M 



Where the 



circular polarization values -f-l and —1. The remaining 
helicity amplitudes may be obtained from parity and per- 
mutation symmetry. Of these five helicity amplitudes, 
three are related by crossing, hence it is sufficient to give 
just three, which are presented in detail in Ref. 

The differential cross section of photon pair production 
from photon fusion, i.e. the box diagram, is 



dcos9 2n s 



(17) 



6 is the scattering angle, a is the fine-structure constant, 
Qf is the charge of the fermion in the loop and the sum 
is over the leptons e and /i and the quarks u, d and s, 
which are the relevant particles in the mass scale that 
we shall discuss. Another possible contribution to this 
continuum process comes from pion loops, which, apart 
from possible double counting, were shown to be negligi- 
ble compared to the above one |2^]. The second sum is 
over the sixteen helicity amplitudes, M\^x^x^\^ , where Ai 
and A2 correspond to polarizations of incoming photons 
and A3 and A4 for the outgoing photons. The matrix ele- 
ments summed over the final polarizations and averaged 
over the initial polarizations is given by 



^|A/p = i{|A/++++P + |Af++__r 
P+4|Af+++_p}. 



2 ' 

|A/,__, 



We consider the scattering of light by light, that is, 
the reaction 77 — > 77, in Au-Au collisions for energies 
available at RHIC, ^/s — 200 GeV/nucleon. We checked 
our numerical code reproducing the many results of the 
literature for the box subprocess, including asymptotic 
expressions for the low and high energies compared to 
the fermion mass present in the loop, and the peak value 
of the cross section (see, for instance, Ref. [pO|). 
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In Fig. (g) the dependence of the ion- ion cross section 
with the cosine of the scattering angle 9, in the two pho- 
ton center-of-mass system, is shown for an invariant pho- 
ton pair mass equal to 500 MeV. It is possible to observe 
that the cross section is strongly peaked in the backward 
direction, but is relatively flat out to cos 9 « 0.4, where it 
starts rising very fast. It is symmetric with respect to 9 
and the same behavior is present in the forward direction. 

It is possible to observe that the cross sections is 
strongly restricted when we introduce an angular cut. We 
will impose in all the calculations throughout this work a 
cut in the scattering angle equal to | cos9\ — 0.5. This cut 
is conservative, but it will make possible to compare the 
cross section of the box diagram with rival processes, that 
will be discussed in the following sections, as well as it is 
enough to eliminate the effect of double bremsstrahlung 
(which dominates the region of | cos9\ w 1). Finally, this 
kind of cut is totally consistent with the requirements 
proposed in Ref. P| 



contribution (up to 0(2 GeV)) are smaller due to their 
masses and charges, because the process is proportional 
to {q'f)'^ where g/ is their charge. Their contribution is 
also insignificant compared to the electron result. 

As discussed in Ref. ^ the 77 scattering can indeed 
be measured in peripheral heavy ion collisions. The cut 
in the angular distribution gives back to back photons in 
the central rapidity region free of the background. How- 
ever, as we shall see in the sequence, there are gaps where 
the 77 77 process is overwhelmed by the presence of 
resonances like the 77, 77/ and others. Even the broad a 
resonance could be of the order of the continuum pro- 
cess. Just to give one idea of the number of events, with 



a luminosity of 2.0 xlO 



[|| and choosing a 
bin of energy of 200 MeV, centered at the energy of 700 
MeV (which is free of any strong resonance decaying into 
two-photons), we have 1532 events/year assuming 100% 
efficiency in the tagging of the ions and photon detection. 
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FIG. 2. Angular distribution of ZZ — > ZZyy scattering at 
an invariant mass of 500 MeV. The scattering angle 9 is in 
the photon pair center-of-mass system. 

The fermions that contribute in the box diagram are 
the leptons e and /i and the quarks u, d and s, which 
are important for the mass range that we are interested 
(heavier quarks will give insignificant contributions and 
the same is true for the charged weak bosons). We as- 
sumed for their masses the following values: TOg = 0.5109 
MeV, ^ 105.6584 MeV, m„ = 5 MeV, = 9 MeV 
and rus = 170 MeV. 

The electron gives the major contribution to the to- 
tal result. The second most important contribution is 
due to the muon, but it is at least one order of magni- 
tude smaller than the electron one. The d and s quark 



B. The process 77 — > i? 77: glueballs 

Photon pair production via the box diagram is a back- 
ground to 77 ^ i? — > 77 process (or vice versa), both 
have the same initial and final states, and for this reason 
they can interfere one in another. Normally the inter- 
ference between a resonance and a continuum process is 
unimportant, because on resonance the two are out of 
phase. 

The total cross section for the elementary subprocess 
77 — *■ i? — *■ 77 assuming a Breit-Wigner profile is 



dM 



167r 



dL 



F2(i?^77) 



dM (Af 2 - )2 + 



41 



^2 • 

total 



(18) 



where M is the energy of the photons created by the 
collision of the ions. T{R — > 77) (= F^-y) and Ttotai are 
the partial and total decay width of the resonance with 
mass mji in its rest frame. 

We are going to discuss only J = resonances made 
of quarks as well of gluons. The reaction 77 — > tt^ — > 77 
was already discussed many years ago where it was 
claimed that the interference vanishes. This result was 
criticized by De ToUis and Violini [Q, affirming (cor- 
rectly) that the interference exists. However, as we will 
discuss afterwards, off resonance the interference is negli- 
gible. If the interference is neg lected, Eq.(|l|) can be used 
and we show in Fig. (||) the result for some resonance pro- 
duction (77, 77', 77(1440), /o(1710)), whose invariant mass 
of the produced photon pair is between 500 MeV and 
2000 MeV. For comparison we also show the curve of the 
continuum process. It is possible to see in that figure 
the well pronounced peaks of the resonances 77 and 77'. 
We assumed for their masses the values of 547.3 MeV 
and 957.78 MeV, respectively, the 77 total decay width is 
equal to 1.18 keV and the 77' one is equal to 0.203 MeV. 
Their partial decay width into photons are 0.46 keV (77) 
and 4.06 keV (77'). 
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The amplitudes for 77 77 intermediated by a pseu- 
doscalar hadronic resonance are : 




600 700 800 900 1000 
M (MeV) 

FIG. 3. Invariant mass distribution of photon production 
(with the cut | cos 91 < 0.5). The solid curve is for the box dia- 
gram, the dashed curves are due to the process 77 — > _R — > 77, 
where R are the pseudoscalars resonances rj and 77' and the 
glueballs candidates ??(1410) and /o(1710). 

We restrict the analysis to the J = glueballs candi- 
dates ?7(1440) and /o(1710). For the r;(1440) we used 
the mass and total decay width values of Ref. [ p3[ , 
niR = 1405 MeV, Ttotai = 56 MeV, for the decay width 
into photons we use the value given in Ref. |^^, T^ry = 5.4 
keV. We see in Fig.(^ that the peak for this resonance is 
of the same order of the continuum process. For the other 
glueball candidate, /o(1710), the peak is clearly above the 
background. For this one we assumed the values listed 
in Ref. of mass and total width, mji = 1715 MeV 
^totai — 125 MeV, and for the two-photon decay width 
we adopted the value encountered by the ALEPH CoUab. 



25 , F, 



77 



21.25 keV. In all these cases the resonances 
can be easily studied in peripheral heavy ion collisions. 

Off resonance we can expect a negligible contribution 
for the process 77 — > i? ^ 77 and consequently the same 
for its interference with the continuum process. How- 
ever, it is instructive to present a more detailed argu- 
ment about why the interference can be neglected. In 
order to do so we are obliged to introduce a model to 
calculate the amplitudes for the process 77 — > i? — > 77. 
These amplitudes will be computed with the help of an 
effective Lagrangian for the pseudoscalar interaction with 
photons (the scalar case will be discussed in the next sec- 
tion), which is given by gpexf^j^pF^f^F^P^p, where gp is 
the coupling of the photons to the pseudoscalar field ^p, 
£\tii/p is the antisymmetric tensor and F^^ is the electro- 
magnetic field four-tensor. 



M 



+++H 

M+-+- 



27r 
27r 
27r 
0, 

-^{F{\st) + F{Xtt) + F{Xn)}, 



F{Xtt), 
F{Xst), 



(19) 



where a is the fine-structure constant, A 
and 



F{x) = 16a; 



r 

2 77 



total 



{mf/mjif 



(20) 



The presence of the fine-structure constant in Eq. ( pj| ) 
is a consequence of the fact that the amplitudes M in 
these equations will be used in Eq. (^7|) , so it is necessary 
to get the correct dependence of the partial cross section 
with this constant. 

A numerical evaluation of the cross section using 
Eq.(|l9|) (for the same resonances present in Fig.(||)) 
shows a totally negligible effect off resonance in compar- 
ison with the box contribution. On resonance the two 
processes are out of phase and the interference is ab- 
sent. We can now proceed with an argument showing 
that the interference is not important. Let us assume 
that off resonance the processes are in phase, and for a 
moment we forget the t and u channels contribution in 
Eq.(|l9[). The s channel contribution can be written as 



M2/ {s — mp, + iT RrnR), and denoting the continuum con- 
tribution by Ml we can write the following interference 
term 

r[{ReMiReM2 + ImMiImlVh) 



,2 \2 



)2 -I- F^,,,.^ 

(s - m%) + {ReMiImM2 - ImMiReM2)T Rma]. 

We can verify that the term proportional to s — inte- 
grates to zero when integrated in a bin centered at to|j. 
With the second term the situation is different. If ImMi 
or ImM2 (assuming ReMi and ReM2 0) then 
there is a nontrivial interference. However, since we are 
dealing with J = amplitudes, the only nonvanishing 
helicity amplitudes are those in which the initial helici- 
ties and the final helicities are equal. Inspection of the 

and AI^^ amplitudes of the box diagram (in 

the limit mf-^rriR) shows that they are purely real, and 
the same happens with M2 (obtained from the s channel 
contribution of Eq.(19)), resulting in a vanishing inter- 
ference! 

Of course, this analysis is model dependent. In partic- 
ular, at the quark level the coupling gp has to be substi- 
tuted by a triangle diagram, which may have a real as 
well as an imaginary part (see, for instance, Ref. p6t). 
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However, this coupling is real for heavy quarks and its 
imaginary part is quite suppressed if the resonance cou- 
ples mostly to light quarks (to/ <C tojj, what happens for 
the resonances that we are considering in the case of the 
u and d quarks, for the s quark the suppression is not so 
strong in the rj case, but we still have an extra suppres- 
sion due to its electric charge). Finally, the interference 
does appear when we consider the amplitudes with the 
resonance exchanged in the t and u channels, but it is 
easy to see that they are kinematically suppressed and 
also proportional to the small value of the total decay 
width. If the total width is large (as we shall discuss in 
the next subsection) the interference cannot be neglected. 
Although the above argument has to rely on models for 
the low energy hadronic physics, we believe that the di- 
rect comparison between the resonant and the continuum 
processes, as presented in Fig.(^), is fairly representative 
of the actual result. 

In Table IV we show the number of events above back- 
ground for the 77, 77' and /o(1710) which, as shown in 
Fig.(^) are clearly above the box contribution. In the 
case of the /o(1710), as well as in the a meson case to 
be discussed in the next subsection, the decay of the res- 
onance into a pair of neutral pions is present. A pair of 
neutral pions can also be produced in a continuous two 
photon fusion process. The rates for all the reactions dis- 
cussed in this work can be modified if both neutral pions, 
no matter if they come from a resonance or a continuous 
process, are misidentified with photons. This accidental 
background can be easily isolated measuring its invariant 
mass distribution, and making a cut that discriminates 
a single photon coming from the processes that we are 
studying, from one that produced two neutral mesons 
(mostly pions) subsequently decaying into two photons. 
For example, in the case of the sigma meson (see next 
section) each one of the neutral pions from its much 
large hadronic decay should be misidentified. These pi- 
ons would produce pairs of photons with a large opening 
angle (j), where cos (0/2) = ^1 — 4to^/to2 . However, 
the calorimeters already in use in many experiments are 
able to distinguish between these two and single photon 
events with high efficiency (see, for instance, Rcf. pTj). 



particle 


events/year 


V 


7.44 X 10'' 


v' 


2.67 X lO"* 


/o(1710) 


42 



TABLE IV. Number of events/year above background for 
the rj, rj' and /o(1710) resonances. 



There is also another accidental background which 
may cause problems to the measurement of the reso- 
nance decay into photons. This is the contribution from 
7— Pomeron^ V 77^. The vector mesons, V, are 
produced with a pt distribution similar to the resonance 
production and at higher rates than some of the processes 
77 —> R, e.g. 7— Pomeron— > ui rate is 10 Hz at RHIC 
pSj , three orders of magnitude higher than for a similar 
mass meson of spin or 2. The lu branching ratio to 
three photons is 8.5%. If a small px photon from this 
decay is undetected, one is left with a low px two-photon 
final state that could be taken for a lighter resonance. At 
higher masses, one also has rjj, 7r°7, K^Ks jX 
as well as 77 /2(1270) tt^tt" and possibly copious 
production of p(1450) and p(1700) by 7— Pomeron inter- 
actions. Clearly a full simulation of all these background 
processes should be kept in mind when measuring two- 
photon final states. 

C. The process 77 ^ ^ 77: the a meson 

The possible existence of light scalar mesons (with 
masses less than about 1 GeV) has been a controversial 
subject for roughly forty years. There are two aspects: 
the extraction of the scalar properties from experiment 
and their underlying quark substructure. Because the 
J = channels may contain strong competing contribu- 
tions, such resonances may not necessarily dominate their 
amplitudes and could be hard to "observe" . In such an 
instance their verification would be linked to the model 
used to describe them. 

Part of the motivation to study the two-photon final 
states in peripheral heavy-ion collisions was exactly to 
verify if we can observe such scalar mesons in its 77 de- 
cay. Although this decay mode is quite rare, it has the 
advantage of not being contaminated by the strong inter- 
action of the hadronic final states. In particular, it may 
allow to investigate the possible existence of the sigma 
meson. This meson is expected to have a mass between 
400-1200 MeV and decay width between 300-500 MeV, 
decaying predominantly into two pions. Of course, an- 
other decay channel is into two photons, with the back- 
ground discussed in the beginning of this Section. 

Recently the E791 Collaboration at Fermilab found a 
strong experimental evidence for a light and broad scalar 
resonance, that is, the sigma, in the 7r~7r"'"7r"'" de- 

cay p9| . The resonant amplitudes present in this decay 
were analyzed using the relativistic Breit-Wigner func- 
tion given by 



m? — ttIq + imoT(m) ' 

with 

^ \PoJ P [Po 
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where m is the invariant mass of the two photons form- 
ing a spin-J resonance. The functions F are the Blatt- 
Weisskopf damping factors [Q: = 1 for spin 

particles, ^F — \l \J\ + (rp*)^ for spin 1 and "^F ~ 
+ 3(rp*)2 + (rp*)''^ for spin 2. The parameter r is 
the radius of the resonance (« 3 fm) |Q and p* = p* (M) 
the momentum of decay particles at mass M, measured in 
the resonance rest frame, pg = pg(m/j). The Dalitz-plot 
of the decay can hardly be fitted without a 0++ {a) res- 
onance. The values of mass and total decay width found 
by the collaboration with this procedure are 478^23 ±17 
MeV and Z^^tf^ ±21 MeV, respectively. 

We will discuss if this resonance can be found in 
peripheral heavy-ion collisions through the subprocess 
77 CT 77. It is important to note that all the values 
related to the cr, like mass or partial widths, that can 
be found in the literature are very different and model 
dependent. In particular, we find the result of the E791 
experiment very compelling and among all the possibili- 
ties we will restrict ourselves to their range of mass and 
total decay width, while we vary the partial width into 
two photons. For the tr decay width into two photons we 
assume the values obtained by Pennington and Boglione, 
3.8 ± 1.5 keV and 4.7 ± 1.5 keV |2|, and the value of 
10±6 keV m. 

It has been verified in the case of tt tt scattering that 
the use of a constant total width in the a resonance shape 
is not a good approximation . In our case we will dis- 
cuss the 77 77 process above the two pions threshold 
where the peculiarities of the broad resonance, basically 
due to the a decay into pions, are not so important. Of 
course, another reason to stay above 300 MeV is that 
we are also far from the pion contribution to 77 scatter- 
ing. In any case we also computed the cross section with 
a energy dependent total width r(m) ~ Fg {v* Ipn?"'^'^ , 
which, as shown by Jackson many years ago |35| ], is more 
appropriate for a quite broad resonance. The net effect is 
a slight distortion of the cross section shape with a small 
increase of the total cross section. Since this one is a 
negligible effect compared to the one that we will present 
in the sequence we do not shall consider it again. 

Off resonance we can expect a negligible contribution 
for the process 77 ^ i? ^ 77 and consequently the same 
for its interference with the continuum process. This is 
true if the resonance has a small total decay width , 
but this is not the u case. To take into account the in- 
terference we must make use of a model to calculate the 
helicity amplitudes of the a meson exchange. Using the 
effective lagrangian ggF'^'''^ F^^^s, where gs is the cou- 
pling of the photons to the scalar field and F'^'' is 
the electromagnetic field tensor the following amplitudes 
comes out ||l|,|2|: 

M++++ = -^F{Xrt), 
M+-+- = -^F{Xtt), 





2-K , , 




M+++ 


- =0, 






- = -'^^{F{\st) + F(\tt) + F(Art)}. 


(21 



where a is the fine-structure constant, A = (m//mj^)^, 
and rt, st and tt are related with the standard Mandel- 
stam variables s, t, and u by r* = -j^, st = and 



These amplitudes and the ones describing the 
fermion (with mass m^) box diagram enter in the 
expression for the differential cross section of pho- 
ton pair production from photon fusion da /d cos 9 = 
(l/27r)(a'*/s)(X;/ q)Y E I^P. to give the total cross sec- 
tion of Fig.(^). We verified that the interference is de- 
structive. 
-4 

10 E . . . . . . 3 
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10 



-9 

10 I • • • • • • ^ 

400 500 600 700 800 900 1000 

U(MeV) 

FIG. 4. Invariant mass distribution of photon pair produc- 
tion. The solid curve is due to box diagram only, the dashed 
one is due to the process 77 —> cr ^ 77 in the Breit-Wigner 
approximation, the dash-dotted is the scalar contribution of 
Eq.(|2^). In all cases F.^^ — 4.7 keV, and the angular cut is 
equal to —0.5 < cos 61 < 0.5 . 

The effective Lagrangian model used to compute the a 
contribution to the photon pair production gives a larger 
cross section than the calculation with the Breit-Wigner 
approximation at energies above M « 600 MeV. It is 
dominated by the s channel contribution. We consider 
the Breit-Wigner result as the best signal representation 
for the resonant process because we are using the E791 
data and this one was fitted by a Breit-Wigner profile. 
The effective Lagrangian gives a nonunitary amplitude 
that overestimates the sigma production above 600 MeV 
and shows the model dependence in the a analysis that 
we commented before. The Breit-Wigner profile is not 
a bad approximation as long as we stay above the two 
pions threshold and in the following we assume that the 
signal is giving by it (the dashed curve of Fig.(^)) and 
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the background is giving by the box diagram result (the 
sohd curve of Fig.(Q)). Note that, due to the destructive 
interference, the actual measurement will give a curve 
below the solid curve of Fig.(||). 

From the experimental point of view we would say that 
the reaction 77 — > 77 has to be observed and any de- 
viation from the continuum process must be carefully 
modeled until a final understanding comes out, with the 
advantage that the final state is not strongly interacting. 
Note that in this modeling the ry meson will contribute 
to 77 ~> 77 in a small region of momentum p9[ | , even so 
it has to be subtracted in order to extract the complete 
a signal. 
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FIG. 5. Significance as a function of decay width into 
two photons, r^.^, for a sigma meson with mass equal to 
478 MeV and total decay width of 324 MeV. The solid 
curve was obtained integrating the cross sections in the in- 
terval 438 < M < 519 MeV, the dashed one in the interval 
300 < Af < 800 MeV. The signal and background are giving 
by the dashed (Breit-Wigner) and solid (box diagram) curves 
of Fig. (^) respectively. 

We changed the values of the a mass and total width 
around the central ones reported by the E791 Collabo- 
ration. We do not observed large variations in our re- 
sult, but noticed that it is quite sensitive to variations 
of the partial decay width into photons. It is interesting 
to look at the values of the significance which is written 
as Ca signal I V Caiiack and characterizes the statistical de- 
viation of the number of the observed events from the 
predicted background. The significance as a function of 
the two photons decay width of the sigma meson, with 
mass equal to 478 MeV and total decay width of 324 
MeV, is shown in Fig.(p|), were we used a luminosity of 
£ = 2.0 X 10^6 cm-2s-~at RHIC and assumed one year 
of operation. The significance is above 2a 95% confidence 
level limit for two photon decay width greater than 4.7 
keV, while for a Scr discovery criteria can be obtained 
with F^-y > 7.5 keV. The numbers in Fig.(|^) were com- 
puted with the signal given by the Breit-Wigner profile. 



and the background by the pure box diagram. The solid 
curve was obtained integrating the cross sections in the 
range of experimental mass uncertainty 438 < M < 519 
MeV, while the dashed curve resulted from the integra- 
tion in the interval 300 < M < 800 MeV. Note that 
there is no reason, a priori, to restrict the measurement 
to a small bin of energy. This choice will depend heav- 
ily on the experimental conditions. Therefore, for values 
of already quoted in the literature the sigma meson 
has a chance to be seen in its two photon decay mode. 
The discovery limits discussed above refer only to a sta- 
tistical evaluation. Our work shows the importance of 
the complete simulation of the signal and background in- 
cluding an analysis of possible systematic errors that may 
decrease the significance. 

IV. CONCLUSIONS 

Peripheral collisions at relativistic heavy ion collid- 
ers provide an arena for interesting studies of hadronic 
physics. One of the possibilities is the observation of light 
hadronic resonances, which will appear quite similarly to 
the two-photon hadronic physics at e"'"e~ machines with 
the advantage of a huge photon luminosity peaked at 
small energies [|l| . Due to this large photon luminosity it 
will become possible to discover resonances that couple 
very weakly to the photons [Q. 

The double Pomeron exchange may be a background 
for the two-photon processes. We have shown that in 
general the Pomeron contribution docs not compete with 
the 77 one, in the case of very heavy ions, after imposing 
the condition for a ultra-peripheral collision. However 
it must be stressed that this comparison has to be per- 
formed for each specific final state, due to the different 
Pomeron coupling to several particles. If the collision in- 
volves heavy ions and the final state we are looking at 
does not have a large coupling to the Pomeron we can 
affirm that double Pomeron exchange can be neglected 
in the evaluation of production rates. In the case of light 
ions Pomeron processes are competitive with the two- 
photons ones, and dominate the cross sections for very 
light ions. 

We discussed the peripheral reaction AA AAjj. 
This process is important because it may allow for the 
first time the observation of the continuous subprocess 
77 — > 77 in a complete collider physics environment. 
This possibility only arises due to enormous amount of 
photons carried by the ions at the RHIC energies. 

The continuous subprocess is described by a fermionic 
box diagram calculated many years ago. We computed 
the peripheral heavy ion production of a pair of photons, 
verifying which are the most important contributions to 
the loop, which turned out to be the electron at the en- 
ergies that we are working, and established cuts that not 
only ensure that the process is peripheral as well as elim- 
inate most of the background. After the cut is imposed 
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we still have thousands of photons pairs assuming 100% 
efficiency in tagging the ions and detecting the photons. 

The continuous 77 77 subprocess has an interest- 
ing interplay with the one resulting from the exchange of 
a resonance. We discuss the resonance production and 
decay into a photon pair. This is a nice interaction to 
observe because it involves only the electromagnetic cou- 
plings of the resonance. Therefore, we may say that it is a 
clean signal of resonances made of quarks (or gluons) and 
its measurement is important because it complements the 
information obtained through the observation of purely 
hadronic decays. It may also unravel the possible amount 
of mixing in some glueball candidates [||. We discuss 
the interference between these process and compute the 
number of events for some specific cases. 

The possibility of observing resonances that couple 
weakly to the photons is exemplified with the a meson 
case. This meson, whose existence has been for many 
years contradictory, gives a small signal in the reaction 
77 — > CT — > 77. However its effects may be seen after one 
year of data acquisition, providing some clue about this 
elusive resonance. Using values of mass, total and par- 
tial widths currently assumed in the literature, we com- 
pute the full cross section within a specific model and 
discuss the significance of the events. Our work shows 
the importance of the complete simulation of the signal 
and background of these processes including an analysis 
of possible systematic errors, indicating that two photon 
final states in peripheral collisions can be observed and 
may provide a large amount of information about the 
electromagnetic coupling of hadrons. 
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